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^ Abstract 

Q^ We discuss the implementation of the "minimal" type III seesaw model, i.e. with one 

^ fermionic triplet, in FeynRules/MadGraph. This is the first step in order to realize a real 

i__i study of LHC data recorded in the LHC detectors. With this goal in mind, we comment on 

^^..^ the possibility of discovering this kind of new physics at the LHC running at 7 TeV with a 

<^ luminosity of few fb~^. 

m 

rn 

O 



> 

X 



^biggio@ifae.es 

^florian. bonnet@pd.infn.it 



1 Introduction 

In a period in which LHC is running and ready to discover new physics, it is of crucial importance 
to have the possibihty of simulating the signals that a particular kind of new physics could give 
in the two main detectors, ATLAS and CMS. In this paper we describe the implementation in 
FeynRules/MadGraph [H [2] of a simple extension of the standard model (SM), the "minimal" type 
III seesaw. This is a first necessary step before performing the analysis of real data, which is the 
ultimate goal of our work and which will be discussed in a future publication. 

As it is well known, oscillation experiments have proved that neutrinos oscillate and therefore are 
massive. However, from the theoretical point of view, the origin of this mass is still unknown. An 
appealing possibilty, also accounting for the smallness of this mass, is the seesaw mechanism: new 
heavy states having a Yukawa interaction with the lepton and the Higgs doublets generate a small 
Majorana mass for the neutrinos, generically suppressed, with respect to charged fermion masses, by 
a factor f /M, where v is the Higgs vev and M the mass of the heavy particle. Depending on the 
nature of the heavy state, seesaw models are called type I [3j , type II |1] or type HI , corresponding 
to heavy fermionic singlet, scalar triplet or fermionic triplet, respectively. If one requires 0{1) Yukawa 
couplings, M should be of the order of the grand unification scale in order to account for neutrino 
masses smaller than the eV. However, in principle the scale can be as low as hundreds of GeV, in 
which case either the Yukawas are smaller or an alternative method, such as for instance an inverse 
seesaw j6] should be at work. In this case the heavy field responsible for neutrino masses could be 
discovered at the LHC. 

As regards collider physics, the seesaws of type II and HI are more exciting, since they can 
be produced via gauge interactions: at difference with singlets, whose production is drastically 
suppressed if the Yukawa couplings are small, triplets can be produced and observed at the LHC if 
their mass is sufficiently small, independently of the size of the Yukawa couplings or mixing angles. 

In the present paper we focus on the type HI seesaw, i.e. the one mediated by fermionic triplets. 
To simplify the implementation of the model in FeynRules, we consider a simple extension of the 
SM obtained by adding a single triplet. Indeed we can safely assume that, unless in case of extreme 
degeneracy, the lightest triplet will be the one most copiously produced and the one which will 
be eventually firstly discovered. In the literature few papers [3 |B1 E] discussing the possibility of 
discovering the type III seesaw at the LHC (at 14 TeV) are present. However so far no code is publicly 
available to perform calculations and simulations in this model. With this paper and the publication 
of the implemented model at the URL http://feynrules.phys.ucl.ac.be/wiki/TypeIIISeeSaw 
we are going to fill this gap. Moreover we briefiy discuss the physics case for LHC running at 7 TeV, 
suggesting that with few fb~^ of luminosity a discovery is already possible. 

This paper is organized as follows. In Sect. 2 the model with the complete Lagrangian and all the 
couplings is reviewed, both in the general and in the simplified case. In Sect. 3 the implementation 
of the model in FeynRules and the checks performed for its validation are discussed. In Sect. 4 the 
physics case at 7 TeV is discussed and in Sect. 5 we conclude. 



2 The model 

The model considered here is the one presented in Ref . [10] . It consists in the addition to the standard 
model of SU(2) triplets of fermions with zero hypercharge, S. In this model at least two such triplets 
are necessary in order to have two non-vanishing neutrino masses. The beyond the standard model 
interactions are described by the following lagrangian (with implicit flavour summation): 

L = TrlEipE] - -Tr[^M^J:^ + !>M*S] - 4)^^V2Y^L - Iv^Fs^S^ , (1) 

with L = (z/, If, (f) = (0+, (f)Y = (0+, {v + H + ir])/y/2f, = zr20*, S^ = CS^ and with, for each 
fermionic triplet, 

D - d -^V2a( ^-'/^ ^^ ] (2) 

Without loss of generality, we can assume that we start from the basis where Ms is real and diagonal, 
as well as the charged lepton Yukawa coupling, not explicitly written above. In order to consider 
the mixing of the triplets with the charged leptons, it is convenient to express the four degrees of 
freedom of each charged triplet in terms of a single Dirac spinor: 

vl/^S+^ + S^. (3) 

The neutral fermionic triplet components on the other hand can be left in two-component notation, 
since they have only two degrees of freedom and mix with neutrinos, which are also described by 
two-component fields. This leads to the Lagrangian 

+ g (w+^^,Pr^ + W^;sf 7^PlvI/ + h.c.) - g W^^^,^ 

- ((j)'f^Yj,UL + y2</)°^ls/L + (jy'-^Y^lL - V2<j)+17^Yi^ + h.c.) . (4) 

The mass matrices of the charged and the neutral sectors need to be diagonalized as they possess 
off-diagonal terms. Following the diagonalization procedure described in Ref. [10], we obtain the 
following Lagrangian in the mass basis: 

^ = ^Kin + ^CC + ^NC + ^NC + ^H + ^H + -^^ + -CJj + ^4,- , (5) 

where 

Ccc = j^{l ^) I'W; [P,gl^ + Pn9TV2) ( ^ ) + h.c. (6) 
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Here Upmns is the lowest order leptonic mixing matrix which is unitary, mi is a diagonal matrix 
whose elements are the masses of the charged leptons, v = a/2(0°) = 246 GeV, e 



'-^Y^M^^Yy, 



^M^^Y^Y^M^^smd 6 



^. The above expressions are all valid at C ( e, e', 6, \/e5, y/e^) 



2.1 The simplified model 

In the previous section the Lagrangian of the type III seesaw model, with a generic number of triplets, 
has been introduced. Since we are interested in LHC physics, we can safely restrict ourselves to the 
case of only one triplet. Indeed, in the presence of more triplets, it will be the lightest the one that 
will be more easily discovered. This will simplify the implementation of the model in FeynRules [^ 
Under this assumption, the new Yukawa couplings matrix reduces to a 1 x 3 vector: 



^s - ( ^s, yj:^ y^-r ) 



(29) 



and the mass matrix M^ is now a scalar. 

The second assumption we will made in the rest of this paper is to take all the parameters real, 
i.e. we do not take into account the phases of the Yukawa couplings nor the ones of the PMNS 
matrix. Barring cancellations, they should not play a role in the discovery process. 

As a consequence e is a 3 x 3 matrix whose elements are 



ea/3 = yM-^ys.>s, , 



(30) 



and e' is now a scalar: 



"'^Notice that while such a simphficd model is appropriate for studies at collider, it accounts only for one neutrino 
mass and therefore does not reproduces the experimental results on neutrino masses. This model should be completed 
with other heavy fields in order to obtain at least two massive light neutrinos. Then this simplified model should be 
viewed as a "low" -energy limit of a more complete theory with heavier states that decouple. If such a hierarchy in the 
masses of the heavy particles is not realized, i.e. if, for example, two or more triplets are degenerate, then the analisis 
will be different. The production cross section for each of the triplet will be the current one, but decays would be 
different, due to the larger number of possibilities for the couplings. 



^ — 2 / ia2 , Az2 



^M^'(Yl+Yl + Yl). (31) 

Finally, we express all the couplings in terms of the mixing parameters, Va = -%M^^Yy;^, since 
they are the parameters which are truly constrained by the electroweak precision tests and the lepton 
flavour violating processes. Then e' = V ■ V'^ while e = V'^ A V . 



By applying these simplifications and redefinitions, the couplings of Eqs. (14)-(28) in terms of 
Ms and Va are obtained; they are shown in Appendix |X] 

3 Implementation of the model in FeynRules and validation 

As discussed in the previous section, the presence of an additional fermionic triplet induces a mixing 
between these new heavy fermions and the light standard model leptons. Then, not only the new 
couplings must be added to the SM Lagrangian, but also SM couplings get modified. In order to 
implement this model in FeynRules, we start from the already implemented SM, contained in the 
file sm.fr, and we add the new coplings and modify the existing ones. The file containing this model 
is named typelllseesaw . f r. In the following we will describe the main features of the implemented 
model, before reviewing the validation checks. 

As shown before, the fermionic triplet can be expressed as a new charged Dirac lepton ^ and 
a Majorana neutral lepton TP . Hence, these two new heavy particles can be viewed as a fourth 
generation in the lepton sector, as suggested by the Lagrangian and couplings written in the previous 
section. Therefore, a new generation index is defined for leptons: 

IndexRange [ Index [LeptonGeneration] ] = Range [4] , (32) 

and charged lepton and neutrino classes have to be extended to include these new heavy parti- 
cles. As for neutrinos, the whole class has to be modified since we are now dealing with Majo- 
rana particles, while in sm.fr the light neutrinos are of Dirac type Q Consequently, the option 

^Note that in the massless hmit the two cases are equivalent. 



Self Con jugate -> True, is turned on. The neutrino class then reads [j 

F[l] == { (33) 

ClassName -> vl, 
ClassMembers -> {vl ,v2,v3,tr0}, 
Flavorlndex -> LeptonGeneration, 
Self Conjugate -> True, 
Indices -> {index [LeptonGeneration] } , 

Mass -> {Mv, {Mvl, 0}, {Mv2, 0}, {Mv3, 0}, {MtrO, 100.8}}, 
Width -> {0, 0, 0, {WtrO, 0.1}}, 

PropagatorLabel -> {"v", "vl", "v2", "v3","tr0"} , 
PropagatorType -> S, 
PropagatorArrow -> Forward, 
PDG -> {8000012,8000014,8000016,8000018}, 
FullName -> {"nul", "nu2", "nu3", "SigmaO"} }. 

Notice that, since neutrinos are Majorana particles, the kinetic term is defined as 

1/2 vlbar.Ga[mu] .del[vl, mu] . (34) 

Analogously, the charged leptons class now reads: 

F[2] == { (35) 

ClassName -> 1, 

ClassMembers -> {e, m, tt,trni}, 
Flavorlndex -> LeptonGeneration, 
Self Conjugate -> False, 
Indices -> {Index [LeptonGeneration] }, 

Mass -> {Ml, {Me, 5.11 * 10(-4)}, {MM, 0.10566}, {MTA, 1.777}, {Mtrch, 101}}, 
Width -> {0, 0, {Wtau, 0.1}, {Wtrch, 0.1}}, 
QuantumNumbers -> {Q -> -1}, 

PropagatorLabel -> {"1", "e", "m" , "tt", "tr-"}, 
PropagatorType -> Straight, 
ParticleName ->{"e-", "m-", "tt-", "tr-"}, 
AntiParticleName -> {"e+", "m+" , "tt+", "tr+"}, 
PropagatorArrow -> Forward, 
PDG -> {11, 13, 15,8000020}, 
FullName -> {"Electron", "Muon" , "Tau" , "Sigma-"} }. 



^The numbers associated to Mass and Width (for S") are variables. 



Notice that the usual PDG codes for hght neutrinos (12, 14, 16) have been replaced by new codes 
(8000012, 8000014, 8000016), since in our model light neutrinos are no longer Dirac particles but 
Majorana ones. Moreover new codes have been provided for the neutral component (8000018) and 
the charged component (8000020) of the triplet. These codes are currently not officially used for 
other particles species and any change should be done very carefully not to interfere with existing 
assignments (see Particle Data Group numbering Scheme pT]). 

Having (re) defined the lepton fields, the interactions can be implemented in the Lagrangian. 
Since the hght leptons couplings to the gauge bosons and Higgs fields are different from the SM 
case, they have been erased and replaced by the ones defined in the previous sections. The matrices 
5'l/R' 9lFw 9l/w 3l/r ^"^^ 9l/r defining the couplings have been introduced as internal parameters 
in order to write the Lagrangian in a clear way. The external parameters, or inputs, are listed in 
Table[lj In this table some values for the parameters of the model implemented in typelllseesaw . f r 
are given, but these are variables that can be modified according to the details of the considered 
model. 

Following the features of the SM implementation, our model presents the characteristic of allowing 
a differentiation between the kinematic mass (or pole mass) of the triplet and the masses entering 
into the couplings definition (equivalent of Yukawa masses). The former are defined under the block 
MASS while the latter are defined under the block NEWMASSES. In particular, for the charged 
fermion masses, we have made the same assignments as in sm.fr: the Yukawa masses for e, fi, u, d, 
s are zero while their pole masses, which are used for example by PYTHIA, are non-zero. This implies 
that any coupling defined in terms of the Yukawa masses will be zero in our model. We have checked 
that turning on this Yukawa masses would amount to a negligible correction. 

3.1 Validation 

In this section we discuss the checks we have performed in order to validate the model we have 
implemented by comparing some numerical results on branching ratios and cross sections obtained 
with typelllseesaw.fr and sm.fr. Moreover, when possible, we will compare the numerical results 
with some analitic expressions. In Table [T]the list of the parameters used for the comparison is given. 

We start by comparing some branching ratios that should not be affected (or very slightly) by the 
presence of the triplet between the FeynRules unitary-gauge implementations in MadGraph/MadEvent 
of the Type III seesaw (typeIIIseesaw_MG) and the SM (siii_FR). These branching ratios have been 
calculated with the program BRIDGE [12j F] and are gathered in Table E^ in Appendix p] They agree 
within 1.5% which roughly corresponds to the intrinsic error of this program; the deviation induced 
by the presence of the triplet is indeed much smaller (~ 0.3%). 

Additionally, these branching ratios can be confronted with the analitic expressions that can be 



''Some care has to be taken when calculating branching ratios with BRIDGE with Majorana particles. Here the 
branching ratios for Z going into Majorana particle has been fixed "by hand" . 



Parameter 


Symbol 


Value in sm.fr 


Value in typelllseesaw 


fr 


Inverse of the electromagnetic coupling 


o^EwiMz) 


127.9 


127.9 




Strong coupling 


as{Mz) 


0.118 


0.118 




Fermi Constant 


Gf 


1.16639e-5 GeV^^ 


1.16639e-5 GeV'^ 




Z pole mass 


Mz 


91.188 GeV 


91.188 GeV 




c quark mass 


rric 


1.42 GeV 


1.42 GeV 




b quark mass 


nib 


4.7 GeV 


4.7 GeV 




t quark mass 


mt 


174.3 GeV 


174.3 GeV 




r lepton mass 


m^ 


1.777 GeV 


1.777 GeV 




Higgs mass 


Mh 


120 GeV 


120 GeV 




Cabibbo angle 


0c 


0.227736 


0.227736 




Electron mass 


rrif. 










Muon mass 


m^ 










Charged heavy fermion mass 


Ms 


- 


101 GeV 




Neutral heavy fermion mass 


Mso 


- 


100.8 GeV 




Light neutrino mass 


mi 












m2 












mg 










PMNS mixing angles 9i2 


C'12 


- 


0.6 




6*23 


023 


- 


0.75 




6'l3 


^13 


- 


0.1 




Heavy-light fermion mixing V^ 


Ve 


- 







y. 


V, 


- 


0.063 




Vr 


Vr 


- 








Table 1: Input parameters for sm.fr and typelllseesaw.fr. 
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Figure 1: Branching ratios of the neutral component (left) and charged component (right) of the 
fermionic triplet in the case Ve = Vr = , V^ = 0.063. The dots correspond to numerically evaluated 
values while the lines correspond to the theoretical predictions. Notice that, as expected from 



Eqs. (36)-(41) in the case of one non-zero mixing angle, the result is the same for charged and 



neutral triplet decay. 



Fig. [T] shows the branching ratios of the charged and neutral component of the fermionic triplet in 
the case Ve = Vr = , V^ = 0.063, while Fig. |2] shows the branching ratios in the case Vr = , Ve = 
V^ = 4.1 ■ 10~^. In both figures, the dots represent the values calculated by BRIDGE while the lines 
correspond to the theoretical predictions. A great agreement is evident. 

Notice that, in case of small mixing angles, the three-body decays of S"*" into E*^ e^(/i^) z/ and 
especially into S'' vr"*" could become relevant jS] and should be taken into account when computing 
branching ratios. We have checked that, for mixing angles of the order of 10^^, Br(S''" — )■ S° vr^) ~ 
10~^, i.e. 2 orders of magnitude smaller than other dominant decays. 

As a second step of the validation procedure, we have computed the cross sections of a selection of 
2 — !■ 2 processes that should not be influenced by the presence of triplets using MadGraph/MadEvent 
and we have compared the results obtained with typeIIIseesaw_MG and sm_FR. Results are gathered 
in Table [S] in Appendix [B] an agreement at the level of 1% is found. 

Finally, we have checked that the production of a pair of triplets at the LHC with a center- 
of-mass energy of 14 TeV obtained with MadGraph/MadEvent matches the previous results in the 
literature [318], see Table [6] in Appendix [b} 
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Figure 2: Branching ratios of the neutral component (left) and charged component (right) of the 
fermionic triplet in the case V^ = , Vg = V^ = 4.1 ■ 10~^. The dots correspond to numerically 
evaluated values while the lines correspond to the theoretical predictions while the lines correspond 
to the theoretical predictions. When both channel with e and fi are open, only one is displayed, 
since, for this particular choice of the mixing angles, they are overlapped. 



4 The minimal type III seesaw model at the LHC at 7 TeV 

4.1 Bounds on the mixing angles 



In Refs. [ini [131 ^M ^^^ bounds on the parameters of the type III seesaw model have been derived. 
The bounds apply to the following combination of parameters: 



— \Y^M-^Y\ „ = \Vo,Vb 

2 I la/3 I " P 



We have then the following constraints: 
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(42) 



(43) 
(44) 
(45) 
(46) 
(47) 
(48) 



Notice that if only Ve or V^ is present the stronger constrain of Eq. (46) does not apply and 
O{10~^) mixings are allowed. On the other side, if both are different from zero, then either one 
of the two is much smaller than the other, effectively reducing this case to the one with only one 
non-zero V^, or they are both O{10~^), in order to satisfy the strong bound of Eq. (46). However, 



as we will discuss later, since the production of the triplet happens via gauge interactions, reducing 



10 



the mixing angle will not reduce the total cross section, so that these bounds have to be taken into 
account, but the mixing angles are not as crucial as in the type I seesaw. 

In this paper we are going to focus on a specific case, in order to illustrate how our model works 
and to show that even with the LHC running at 7 TeV there is the possibility of testing the low 
scale type III seesaw. We are going to give the cross section of the relevant channels for the case 
Ve = Vr = , Vfj_ = 0.063. This case corresponds to the maximum allowed mixing angles. If the 
mixing is so large, then some cancellation or an extended seesaw mechanism like the inverse seesaw 
must be invoked in order to obtain the correct value for neutrino masses. However, all the discussion 
we perform in this section applies also in the case of small mixing. In the next sections we are going 
to discuss the triplet production and decays, give the cross sections which are relevant for discovery 
and discuss the main backgrounds which affect the measurement and the main cuts that could be 
implemented in order to reduce it. A more detailed study is beyond the scope of this work. 

4.2 Triplet production and decay 

At the LHC triplets are mainly produced in pair. In Table |2] production cross sections for different 
mass values are collected, with the acceptance cuts listed in Table |3j Since the triplets are produced 
via gauge interactions, the production cross sections do not depend on the mixing parameters. After 



production, the triplets decay inside the detector according to the expressions displayed in Eqs. (36) 



(41). While the decay width depends strongly on the value of the mixing angles Va, the branching 
ratios dependence is very mild. Since we are always in the narrow width regime, the total cross 
section is driven only by the mass of the triplet (for the production) and its branching ratios (for 
the decays). Therefore, a non-discovery at the LHC will permit to constrain the mass of the triplet, 
after some assumption on the branching ratios have been done. 

Once the triplets have decayed into leptons and gauge bosons, the latter will then decay into 
charged leptons, quarks, which will show up as jets (and leptons, when heavy quarks decay semilep- 
tonically), and neutrinos, which will manifest themselves as missing energy. Final states can be 
classified according to the number of charged leptons. The type HI seesaw can give rise to final 
states with up to 6 leptons. However, it has been shown that the cross sections for 6-, 5- and 4- 
leptons final states is to low for being useful for discovery, already at 14 GeV 0; therefore, we will 
not consider them here [_| On the other hand, the most promising channels are the 3-leptons and 
the dileptons, i.e. with 2 leptons of the same sign. In the following sections we are going to discuss 
these channels and the main backgrounds which affect them. 

^However, since the probability of missing a lepton is relatively high for multilepton channels, when generating 
events to study the possibility of having a signal in the 3- and 2-leptons channels, events with 4 leptons should be 
generated too. The inclusive 4-leptons final state cross section varies between 10-20 fb for triplet masses in the range 
100-140 GeV. 



11 



Ms 


a(pp->S+S0)(/6) 


a(pp^S+S-)(/6) 


a(pp^S-S0)(/6) 


100 


4.329e+3 


3.339e+3 


2.325e+3 


120 


2.157e+3 


1.629e+3 


1.106e+3 


140 


1.200e+3 


8.882e+2 


5.894e+2 


160 


7.215e+2 


5.229e+3 


3.387e+2 


180 


4.555e+2 


3.249e+2 


2.059e+2 


200 


3.006e+2 


2.109e+2 


1.311e+2 


300 


5.488e+l 


3.580e+l 


2.027e+l 


400 


1.434e+l 


8.777 


4.632 


600 


1.527 


8.576e-l 


4.118e-l 


800 


2.097e-l 


1.132e-l 


5.139e-2 


1000 


3.133e-2 


1.774e-2 


7.401e-2 



Table 2: Production cross sections at 7 TeV. 



Acceptance Cuts 


Pt^ > 20 GeV 


7]j<5 


ARjj > 0.001 


PTi > 10 GeV 


Vi < 2.5 


ARu > 



Table 3: Acceptance cuts used for production simulations at 7 TeV and 14 TeV. 

4.3 The most relevant final states 

Tables [7] and |8] in Appendix [C] display the cross sections for the intermediate and final states with 2 
and 3 leptons at different mass energies. |j While the intermediate ones are calculated with MadGraph, 
the final ones are obtained by multiplication with the corresponding branching ratios. From a quick 
look to these tables one can see that even with LHC running at 7 TeV, with the few fb~^ of luminosity 
which are expected to be reached by the end of 2011, several events are expected, for low triplet 
mass. In the 3-leptons table, in the total cross section we have isolated the channels with leptons 
not-coming from Z decay. Indeed, when the cut on the invariant mass of the leptons will be applied 
in order to reduce the background events coming from Z decay (see later), these events will mostly 
disappear. Then the numbers we quote in blue in Table |8] can be considered the effective cross section 
after the application of this cut. 

By looking at these table we see that there are 4 possible final states with 2 and 3 leptons: 

A) 3 leptons + missing transverse energy (MET); 

B) 3 leptons + 2 jets + MET; 

^We give numbers for the case of mixing with muons exclusively, however similar results apply when the final 
states contains electrons as well. On the other hand, they do not apply completely to taus. Indeed, taus are not 
detected as such, because of their fast decay. Moreover, in a detector like CMS, leptons coming from taus decay are 
not distinguished from prompt leptons and therefore identified taus are only hadronic taus. 
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Figure 3: Dominant process for the discovery channel for the fermionic triplet at the LHC 



C) 2 same-sign leptons + 4 jets; 

D) 2 same-sign leptons + 2 jets + MET. 



In what follows we are going to discuss the main features of all of them. We have simulated pp — )■ 
S+S'^ — 7- fi~^ij,'^n~ + usl+jets) with MadGraph/MadEvents, hadronization being obtained with the 
help of PYTHIA [15]. The CMS detector has been simulated via the PGS software [16j. 




50 100 150 200 250 300 350 400 



Figure 4: Invariant mass of the two /i+ for a luminosity of 30fb and Ms = 100 GeV. Pre-selection 
cuts selected only the events with 3 charged leptons among which 2 positive muons. 



leptons -|- MET. This is probably the best discovery channel: indeed the background is more 
easily reduced due to the absence of jets in the final state. The dominant process generating 
it is depicted in Fig. |3| In an ideal detector where jets are not misidentified with leptons, the 
only background sources would be WW, WWW, WZ and ZZ when a lepton is missed. In 
practice jets should be added to these background; however, as it is discussed later, all these 
background should be under control. 
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In this channel, the invariant mass m^+^+ of the two same-sign muons presents a long tail in 
the high energy region that is characteristic of the presence of new physics, see Fig. |4[ and 
can be exploited to reduce the background. Moreover, this is typical of this kind of seesaw, 
permitting thus to distinguish among type I, II and III [7J. 



0.12 



0.08 




0.02 



20 40 60 80 100 120 140 



Figure 5: pt distribution of the different leptons for M^ = 100 GeV. The black, red and blue curves 
represent the lepton with the highest, intermediate and smallest pt respectively. Pre-selection cuts 
selected only the events with 3 charged leptons among which 2 positive muons. 



3 leptons + 2 jets + MET. This channel is probably the best one in order to reconstruct the 
mass of the triplet. Moreover it can be used also to discriminate between type II and type III 
seesaw [7]. It also appears in the type I seesaw with a gauged U{1)b~l [T7]. In this case the 
reduction of the background can be more complicated, due to the impossibility of applying a 
jet veto. Essentially all the sources listed in the next section constitute a background for this 
channel. A precise estimation of the sensitivity to this new physics would require the complete 
simulation of the background and a detailed analysis, which is beyond the scope of this work. 
However, we will show later that the possibility of reducing the background to "reasonable" 
levels is realistic. 

Once the triplet has been observed, its mass needs to be measured. To this aim, this channel, 
emerging from the process pp — )■ (S^ — )■ i."^ Z / H){TP — )■ l^W^) with Z/H decaying into jets, is 
the best one. Indeed the momentum of the Z/H boson is reconstructed from the jets momenta, 
while its combination with the momentum of one of the two same-sign leptons gives the mass 
of the charged triplet. Since there are two possibilities for this combination, the chosen one 
will be that giving closest invariant mass for the recontructed charged and neutral triplets, 
where the latter is given by the combination of the momenta of the two remaining leptons plus 

metQ 



'^The neutrino longitudinal momento should be added as well [7]. 
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The reconstructed mass of the charged and neutral triplet are shown in Fig. |6] where no cuts 
has been applied. Note that a selection cut on the invariant mass rrijj of the jets 



\m 



']] 



M2/^|<10GeV (49) 

will improve the mass reconstruction. Even if the background is added, a clear peak in the 



— lOOGeV 
— 140 GeV 




— lOOGeV 
— 140 GeV 




50 100 



Figure 6: Reconstructed mass of the charged triplet (left) and neutral triplet (right), for a luminosity 
of 30fb~^, in the case M^ = 100 GeV (black curve) and Ms = 140 GeV (red curve). Pre-selection 
cuts selected only the events with 3 charged leptons and at least 2 jets. 



reconstructed mass will still be visible, which should also permit to distinguish from type II 
seesaw [7j. 

same-sign leptons + jets (+MET) As it is clear from Table [TI the cross section for these 
final states are quite large, even larger than the ones for 3 leptons final states. However 
here jets are always present, which can render a bit more difficult the background reduction. 
The backgrounds are essentially the same as in the previous channel and indeed it has been 
shown [7] that the discovery and the discriminatory potentials of the 2- and 3-leptons final 
states are similar too. A realistic study, especially a study on real data, should consider this 
channel as well. 



4.4 Background 

The main background sources for the channels discussed above are : tt, ttW , WW, WZ, ZZ, Zti, 
Zbb and 3 gauge bosons. The same background plus additional jets should be considered as well, 
both if looking at final states with jets or no: some jets can be indeed misidentified as leptons. In the 
following we will give a brief description of each background and of the cuts that can be implemented 
in order to reduce it. Whenever the cross section for the different background under study has not 
been measured, we have used MadGraph/MadEvent to obtain the cross-sections for LHC running 
at 7 TeV and compared our results with previous results obtained by the CMS collaboration 
whenever possible. All backgrounds have been simulated with and 1 additional jets. 
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tt. The production of a pair of top quarks decaying into bW, one of the b giving a lepton and 
the W decaying leptonically, is a source of background with a large cross section. At 7 TeV 
the production of a top quarks pair has been measuread by CMS [IH] and ATLAS [20] to be 
(jjj = (173^32) pb and (171 ± 20 ± 14^g) pb, with an integrated luminosity of 36 and 35 pb~^, 
respectively. Combining the branching ratio BR{W — )■ lu) = 30% with the 10% of branching 
ratio for the semileptonic decay of the b, the final cross section for such background should be 
around 0.15 — 1.5 pb depending on how many different lepton flavors one expect in the flnal 
state. In the case where the signal flnal state does not contain jets (at the parton level), a 
cut on the number of jets will reduce this background to negligible levels. 6-tagging could be 
applied in order to reduce it when channels with jets are considered. 

ttW. Here the two tops decay into a W plus jets. The third W ensures the presence of three leptons 
in the flnal state. The presence of jets makes this background negligible when looking to three 
leptons + MET without jets. On the other hand, when channels with jets are considered, 
this background should be carefully studied. We found auwj ~ 230 fb. The production cross 
section for ttW should then be larger, but considering the appropriate branching fractions, the 
flnal cross sections should be of few fb, depending on the number of jets. 

WW. This is a large source of background. At 7 TeV, it has been measured by CMS [21] and AT- 
LAS [22J tobe: aww = 41.1 ±15.3(stat.)± 5.8(syst.) ±4.5(lumi.) pb and av^vK = 41+^^(stat.)± 
5(syst.) ± l(lumi.) pb, with an integrated luminosity of 36 and 34 pb~^, respectively. CMS 
collaboration also found [23| : a{pp — >■ VTW^ + X) = 55.3±3.3(stat.)±6.9(syst.)±3.3(lumi.) pb. 
But pre-selection cuts (3 charged leptons out of which 2 have the same sign, 2 hard leptons) 
should reduce it to a negligible level. 

WZ. The CMS collaboration measured [23] : a{pp -^ WZ + X) = 17.0 ± 2.4(stat.) ± l.l(syst.) ± 
1.0(lumi.) pb. This will give ~ 60 fb for the flnal state cross section. A cut on the invariant mass 
of two leptons with opposite sign, \Mz — mii\ > 10 GeV , can be applied in order to eliminate 
leptons coming from Z decay. Moreover, if one considers leptons with different flavour, like for 
instance the channel e~ fi~^fi~^ + MET, this will be free from such a background. 

ZZ. This channel is a background when one of the lepton is lost. It has been measured at the LHC 
by the CMS collaboration [23] : ^ipp ^ ZZ + X) = 3.8l}J(stat.) ± 0.2(syst.) ± 0.2(lumi.) pb. 
Again, cuts on the invariant mass of opposite signs leptons should allow to reduce it to a 
negligible level. 

ttZ and bbZ. These constitute a background for flnal states involving jets. The production cross 
section is relatively large: a^jz = 205 fb and a^j^ = 50 pb. However, the cuts on the invariant 
mass of the leptons as well as 6-tagging should reduce them to negligible levels. 

WWW. Among the 3 gauge bosons background, this is the one with highest cross section. The 
production cross-section for three W bosons is anyway lower than other background considered: 
(Jwww = 71 fb, which becomes really negligible when the flnal state is considered. 

All theses background sources can be reduces by cuts on the pt of the leptons which are hard in 
the signal flnal state. Additional cuts on number of jets or opposite-sign leptons' invariance mass 
can further help to improve the signal over background ratio. 
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As it is clear, the aim of this section was just to describe the main backgrounds affecting the 
considered signals. In order to give precise estimation the entire simulation of the background should 
be performed. 

4.5 Other relevant cases 

Even if we have discussed in details only the case of large mixing with muons, there are other cases 
which can be relevant. Here we briefly sketch their characteristics. 

Mixing with electrons or taus. As already discussed in the literature [7j, the situation for mixing 
with electrons is similar to the one with muons and our analisis can be applied to it as well. 
On the other side, since detecting taus is more complicated, the discovery potential of channels 
involving taus is believed to be smaller. 

Mixing with 2 or 3 charged leptons. In such a case the triplet can couple to more than one 
family. The mixing angles are thus more constrained. As we have already shown (see Figs, fl]- 
|2]), the simultaneous presence of two (or three) non zero Va would reduce the corresponding 
branching ratio by a small factor: if, for instance, two of them are taken to be equal, then the 
corresponding branching ratio will be decreased by a factor 2 with respect to the case with 
only one non-zero mixing angle (see Figs. fl](2]). However the pair production cross section of 
triplets is not affected by the mixing values and thus only the branching ratios and the mass 
of the triplet drive the relevant processes studied here. 

Small mixing angles, O{10~^). This case is the "most natural" one, since here small neutrino 
masses can be accomodated without any cancellation or further source of suppression Q Such 
small mixing angles drastically reduce the value of the triplet decay width, so that displaced 
vertexes up to few millimeters can be present (see also [S]). In case of finding an excess of 
events in some of the considered channels, the measurement of these displaced vertexes could 
be a clear signal that we are in presence of this kind of physics. The possible presence of a 
displaced vertex have to be taken into account when defining the reconstruction parameters for 
the data analisis (for example to reconstruct an interaction vertex). A detailed study of this 
topic is postponed to the analisis of real data. A part from this, in general the cross sections 
are not affected and the analisis can proceed as in the case of large mixing. 



5 Conclusions 

In this paper we have described in details the minimal type HI seesaw model and its implementation 
in FeynRules/MadGraph. In particular we have explicitly written all the couplings and we have 
discussed the tests we have performed in order to validate the implemented model. Even if the 



^Notice that in this case the approximation of taking zero neutrino masses is no longer consistent and they should be 
turned on in the numerical simulations; for consistency also non-zero electron and muon masses should be considered, 
even if the effect of all these masses turns out to be negligible. 
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model has been tested only with MadGraph which uses the unitary gauge, the Goldstone bosons have 
been implemented as well, so that it can be used also with other Monte Carlo generators such as 
CalcHep [21]. As already stressed in the Introduction, this is a necessary step to be done before 
proceeding to the analysis of real LHC data. 

In order to show an example of the utility of our model, we have focused on a particular case 
-large mixing with muons, V^ = 0.063, and small triplet masses, 100 GeV, 120 GeV, 140 GeV- and 
for these cases we have calculated the cross sections of the relevant channels at the LHC running 
at 7 TeV. We have shown that several events are expected for a luminosity of few fb~^. We have 
discussed the main background sources and the methods that can be employed in order to reduce it. 
A more detailed study is beyond the scope of this work, but, still at this level, we can expect that a 
discovery at the LHC is possible, even in the 2011 run, if the mass of the triplet is low enough and 
the background rejection is good. Otherwise, in case of non-discovery, an upgrade of the bounds on 
the triplet mass can be set. 
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A The explicit Lagrangian in the minimal model 
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In the above expressions repetead flavour indexes are summed. As we will discuss later, we will 
take neutrino masses equal to zero, except in the case of small mixing angles [^ 



^In this case, indeed, for consistency we will turn neutrino masses, as well as electron and muon masses, on. 
However, this will not basically affect the result. 
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B Tables for the validation of the implementation 



Process 


sm_FR 


typeIIIseesawl_MG 


comparison 


top decay 


1.53174916 


1.55409729 


1.45899% 


W decay 


2.00335798 


2.00322925 


0.00642571% 


Z decay 


2.41539342 


2.41481975 


0.0237506% 


BR(w+ -> V e+ ) 


1.11025062e-01 


1.11142e-01 


0.105326% 


BR(w+ -> V m+ ) 


1.11036355e-01 


1.11331e-01 


0.265359% 


BR(w+ ^ V tt+ ) 


1.12013868e-01 


1.11018e-01 


0.8962% 


BR(w+ ^ c d ) 


1.69615944e-02 


1.69574065e-02 


0.0246905% 


BR(w+ ^ u d ) 


3.14853587e-01 


3.16304871e-01 


0.460939% 


BR(w+ ^ c s ) 


3.17238100e-01 


3.16278512e-01 


0.302482% 


BR(w+ -> u s ) 


1.68714343e-02 


1.69683505e-02 


0.574441% 


BR(z ^ e- e+ ) 


3.45878542e-02 


3.45049797e-02 


0.239606% 


BR(z -)■ m- m+ ) 


3.46182266e-02 


3.49703234e-02 


1.01709% 


BR(z -^ tt- tt+ ) 


3.45433552e-02 


3.45770661e-02 


0.0975901% 


BR(z -)■ invisible) 


0.205237 


0.205557 


0.155917% 


BR(z ^ b b ) 


1.51238258e-01 


1.50200176e-01 


0.686388% 


BR(z ^ c c ) 


1.17361782e-01 


1.17167722e-01 


0.165352% 


BR(z -> d d ) 


1.52782011e-01 


1.52925551e-01 


0.0939509% 


BR(z ^ s s ) 


1.52615959e-01 


1.51787006e-01 


0.543163% 


BR(z ^ u u ) 


1.17015696e-01 


1.18309630e-01 


0.10578% 



Table 4: Comparison of decay widths and branching ratios between the model sm_FR and 
typeIIIseesawl_MG. 
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Process 


sin_FR 


typeIIIseesawl_MG 


comparison 


e+e' — )> e+e" 


7.457e+2 


7.450e+2 


0.095 % 


e + e^ -> ;U + yU" 


1.125e-l 


1.126e-l 


0.09 % 


e+e" — ;> z/+z/" 


5.185e+l 


5.180e+l 


0.10% 


T + T- -^ W + W- 


2.629e+0 


2.625e+0 


0.15% 


T + T- -^ ZZ 


1.448e-l 


1.449e-l 


0.07% 


T + T- -^ Z7 


7.208e-l 


7.219e-l 


0.15% 


T + T^ — 7- 77 


1.020e+0 


1.020e+0 


- 


ZZ ^ ZZ 


5.997e-l 


5.996e-l 


0.017% 


W+W- -^ ZZ 


2.996e+2 


2.995e+2 


0.033% 


HH ^ ZZ 


6.763e+l 


6.763e+l 


- 


HH -^ W+W- 


1.046e+2 


1.039e+2 


0.57% 


GG^GG 


3.084e+5 


3.079e+5 


0.16% 


uu -)> GG 


1.981e+2 


1.980C+2 


0.05% 


uu -> W+W- 


8.711e-l 


8.720e-l 


0.10% 


uu — )■ ZZ 


8.783e-2 


8.800e-2 


0.19% 


uu — )■ ^7 


1.215e-l 


1.216e-l 


0.08% 


MM -> 77 


6.725e-2 


6.714e-2 


0.13% 


UU — )■ ss 


7.809e+0 


7.807e+0 


0.026 % 


■urf — )■ cs 


1.040e-l 


1.040e-l 


- 


MS — )■ cd 


3.000e-4 


2.999e-4 


0.033% 


tt^GG 


7.352e+l 


7.349e+l 


0.027% 


ti -^ w+w- 


7.521e+0 


7.512e+0 


0.12% 


ti^ ZZ 


7.875e-l 


7.899e-l 


0.30% 


tt^ Z-i 


4.778e-l 


4.771e-l 


0.15% 


ti — )■ 77 


3.096e-2 


3.091e-2 


0.161% 


tt -> uu 


3.139e+0 


3.130e+0 


0.28% 



Table 5: Selection of 2 — )> 2 processes. The FeynRules generated Standard Model implemen- 
tations in MadGraph/MadEvent is denoted sin_FR and the one of the type III Seesaw is denoted 
typeIIIseesawl_MG. The center-of-mass energy is fixed to 1 TeV and a pt cut of 20 GeV is applied 
to each final state particle. 
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Ms 


a(pp->S+S0)(/6) 


a(pp^S+S-)(/6) 


a(pp^S-S0)(/6) 


100 


1.126e+4 


9.125e+3 


6.914e+3 


120 


5.818e+3 


4.673e+3 


3.480e+3 


140 


3.373e+3 


2.673e+3 


1.957e+3 


160 


2.100e+3 


1.646e+3 


1.184e+3 


180 


1.382e+3 


1.071e+3 


7.604e+2 


200 


9.471e+2 


7.273e+2 


5.073e+2 


300 


2.136e+2 


1.564e+2 


1.023e+2 


400 


7.012e+l 


4.847e+l 


3.039e+l 


600 


1.280e+l 


8.307 


4.713 


800 


3.290 


1.993 


1.068 


1000 


1.018 


5.896e-l 


2.978e-l 



Table 6: Production cross sections at 14 TeV. These values have been obtained with 
MadGraph/MadEvent and the acceptance cuts implemented are listed in Table |3} Fig. [7] shows the 
interpolated curves. 



CT(pp-ZZ) (fb) 



1000 







1000 
Mz (GeV) 



Figure 7: Production of a pair of triplets at 14 TeV at the LHC The mixing parameters as been 
set to Vf, = 0.063 and K = K = 0. 
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C Cross sections of the relevant channels at 7 TeV 



Process 


Cross Sections (fb) 


Final State 


Final State Cross Section (fb) 




100 GeV 


120 GeV 140 GeV 




100 GeV 


120 GeV 


140 GeV 


Final State ++ 


W-^i+Z^i+ 


2.36e + 2 


2.02e + 2 


1.16e + 2 


fi'^ fi~^ hadr 


108 


92.7 


53.4 










^'^H^i/vhadr 


32.4 


27.8 


15.9 


W^li+W+u 


1.66e + 3 


6.06e + 2 


2.82e + 2 


^'^^'^uuhadr 


124 


45.3 


21.1 


W-^i+h^i+ 


1.22e-3 


1.39e-l 


1.40e + l 


^'^^'^hadr 


- 


- 


8.9 










fi'^fi'^uuhadr 


- 


- 


- 


Total Cross Sections ^'^ ^'^ + jets + missing Et 


156.4 


73.1 


37.0 


Total Cross Sections /i^/i"^ + jets 


108 


92.7 


62.3 


Final State 


w+ii-z^i- 


1.27e + 2 


1.04e + 2 


5.67e + 1 


fi~fi~hadr 


58.3 


47.7 


26.1 










^~ ^~vvhadr 


17.4 


14.3 


7.8 


W+iJ-W-u 


8.94e + 2 


3.11e + 2 


1.39e + 2 


fi~fi~uuhadr 


67.0 


23.3 


10.4 


W+i^-hfi- 


5.87e - 6 


7.13e-2 


6.86 


fi~fi~hadr 


- 


- 


4.4 










^~ ^~vvhadr 


- 


- 


- 


Total Cross Sections ^~ ijT + jets + missing Et 


84.4 


37.6 


18.2 


Total Cross Sections f^~fi~ + jets 


58.3 


47.7 


30.5 



Table 7: Final states with two muons of the same sign for V^ = V^ = 0, V^ = 0.063. The final 
cross sections have been computed using the measured branching ratios, except for the Higgs, whose 
branching ratios have been calculated assuming a mass of 120 GeV. Only channels with a final cross 
section higher than 0.1 have been reported. 
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Process 


Cross Sections (fb) 


Final State 


Final State Cross Section (fb) 




100 GeV 


120 GeV 140 GeV 




100 GeV 


120 GeV 


140 GeV 


Final State + + - 


W+n-W+u 


1.66e + 3 


6.08e + 2 


2.82e + 2 


jX^ H'^ li~ wp 


20.9 


7.7 


3.5 


W'I2+W+U 


1.66e + 3 


6.06e + 2 


2.82e + 2 


jl^jJi^ll^VVV 


20.9 


7.7 


3.5 


W'^fi-Zfi+ 


2.36e + 2 


2.03e + 2 


1.16e + 2 


11^ 11^ li~ vhadr 


18.2 


15.7 


8.9 










li^ ix^ ix~vvv 


5.5 


4.7 


2.7 


W-i^+Z^i+ 


2.36e + 2 


2.02e + 2 


1.16e + 2 


fi'^fi'^H~i'hadr 


18.3 


15.6 


8.9 










fl'^fI~^fi~l'UU 


5.5 


4.6 


2.6 


W+uZu 


4.62e + 2 


4.02e + 2 


2.32e + 2 


11^ 11^ ix~vvv 


1.8 


1.6 


0.9 


Z/i+Zu 


6.55e + 1 


1.35e + 2 


9.48e + 1 


11^ 11^ ix~ vhadr 


1.6 


3.2 


2.3 










ix^ ix^ ix~vvv 


0.47 


0.98 


0.68 


Zfi^hu 


6.80e - 4 


1.54e - 1 


2.28e + l 


fj,'^ fj,'^ fj,~ uhadr 


- 


- 


0.76 


W-uZ^i+ 


3.61e + 2 


3.08e + 2 


1.71e + 2 


^'^ ^'^ fx~ uhadr 


8.4 


7.2 


4.0 


W+n-hiJ+ 


1.22e-3 


1.39e-l 


1.40e + l 


fi^ fi^ fi^ uhadr 


- 


- 


1.5 


W-fi+hfi+ 


1.22e-3 


1.39e-l 


1.40e + 1 


fi^ fi^ fi^ uhadr 


- 


- 


1.5 


Total Cross Sections /i+zi^/i" + jets + missing Et 


46.5 


41.7 


27.9 


Total Cross Sections ji^ fi'^ ^~ + jets + missing Et (only via W) 


36.5 


31.3 


20.8 


Total Cross Sections /i"'"/i "•"//" + missing Et 


55.1 


27.3 


13.9 


Total Cross Sections /i"'"/i"'"/i~ + missing Et (only via W) 


52.8 


24.7 


12.3 


Final State H 


W'li+W-u 


8.96e + 2 


3.13e + 2 


1.39e + 2 


H~H~H~^i'i>i> 


11.2 


3.9 


1.7 


W+fi-W-u 


8.94e + 2 


3.11e + 2 


1.39e + 2 


fl~fl~fl~^h'UU 


11.1 


3.9 


1.7 


W-fi+Zfi- 


1.27e + 2 


1.04e + 2 


5.67e + 1 


fi^fi^fi^uhadr 


9.8 


8.0 


4.4 










^~ ^~ ^^vvv 


2.9 


2.4 


1.3 


W+f^-Zfx- 


1.27e + 2 


1.04e + 2 


5.67e+l 


fi~fi~fi~^i'hadr 


9.8 


8.0 


4.4 










fl~fl~fI~^l'UU 


2.9 


2.4 


1.3 


W-vZv 


2.49e + 2 


2.07e + 2 


1.13e + 2 


fi~fi~fi~^i'i'i' 


1.0 


0.8 


0.4 


Z^i^Zv 


3.53e + 1 


6.93e + 1 


4.65e + 1 


fi^fi^fi'^uhadr 


0.85 


1.7 


1.1 










fi^fi^fi^uuiy 


0.25 


0.5 


0.3 


Zn^hv 


3.27e - 4 


7.87e - 2 


1.12e + l 


fi^fi^fi'^uhadr 


- 


- 


0.37 


W^vZ\r 


3.62e + 2 


3.07e + 2 


1.72e + 2 


fi^fi^^'^uhadr 


8.4 


7.2 


4.0 


W-^i+h^i- 


5.87e - 4 


7.13e-2 


6.86 


fi~fi~^~^i'hadr 


- 


- 


0.7 


W+fi'hfi- 


5.86e - 4 


7.10e-2 


6.87 


fi^fi^fi'^uhadr 


- 


- 


0.7 


Total Cross Sections /i+/x~/i~ + jets + missing Et 


28.9 


24.9 


15.7 


Total Cross Sections fi^fi^ fi~ + jets + missing Et (only via W) 


19.6 


16.0 


10.2 


Total Cross Sections n'^n~n~ + missing Et 


29.4 


13.9 


6.7 


Total Cross Sections fi'^fi~fi~ + missing Et (only via W) 


28.1 


12.6 


6.0 



Table 8: Final states with three muons for Ve = Vr = 0, V^ = 0.063. The final cross sections have 
been computed using the measured branching ratios, except for the Higgs, whose branching ratios 
have been calculated assuming a mass of 120 GeV. Only channels with a final cross section higher 
than 0.1 have been reported. As for the total cross sections, we have isolated the ones where the 
muons are generated via W decay, since almost al|4he muons generated via Z decay will be removed 
by the cut implemented to reduce the Z background. 



